Cross phase modulation in a five— level atomic medium: Semiclassical theory 
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The interaction of a five-level atomic system involving electromagnetically induced transparency 
with four light fields is investigated. Two different light-atom configurations are considered, and 
their efficiency in generating large nonlinear cross-phase shifts compared. The dispersive properties 
of those schemes are analyzed in detail, and the conditions leading to group velocity matching for 
two of the light fields are identified. An analytical treatment based on amplitude equations is used 
in order to obtain approximate solutions for the susceptibilities, which are shown to fit well with 
the numerical solution of the full Bloch equations in a large parameter region. 
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I. INTRODUCTION 

An cfRcient cross-phase modulation (XPM) in quan- 
tum and semiclassical regimes is both interesting and 
useful in many possible applications, such as those in 
optical communications 0, optical Kerr shutters 0, 
quantum non-demolition measurements Q and quantum 
phase gates 0. In all of these, but the last two especially, 
a large XPM is desirable for low pump powers and high 
sensitivities. 

In a standard three-level cascade scheme, shown in 
Fig. njj, nonlinear effects are obtained alongside absorp- 
tion, which increase as the fields are tuned closer to the 
atomic transition Q . To reduce the absorption to an ac- 
ceptable level, light fields need to be strongly detuned 
from the intermediate atomic level |2), simultaneously 
reducing however the size of the nonlinearity, since both 
are inversely proportional to the square of the detuning. 

Extensive studies aimed at avoiding this problem have 
been performed in recent years. A promising candidate 
emerged with the use of quantum coherence effects in 
the interaction of light with multilevel atoms. Coherent 
population trapping (CPT) j^] and in particular the re- 
lated effect of electromagnetically induced transpare ncy 
(EIT) 0, Q have been studied theoretically [1 H Eg 
and experimentally 0, IllL [12. I'i] in various energy-level 
schemes based on a generic A-scheme (see Fig. ^) . At 
resonance (Si = in Fig. ^), the presence of the cou- 
pling field (with Rabi frequency 0,2) cancels, by destruc- 
tive interference, the absorption on the probe transition 
(with Rabi frequency ^1), and renders the medium trans- 
parent for the probe beam. A more general condition for 
EIT is two-photon resonance, a condition that is satisfied 
when the frequency difference between the fields matches 
the energy gap between levels |1) and |3). However, on 
the exact EIT resonance, probe field decouples from the 
atoms, making the dynamics purely linear. 

Optical nonlinearities in a multilevel atomic or molec- 
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FIG. 1: (a) Three-level A-scheme for EIT. Transitions are 
driven by the probe and coupling fields, with Rabi frequency 
f^i and Q2 respectively. When the probe detuning 5i matches 
the two photon Raman-resonance condition with the coupling 
field, the atomic medium becomes transparent for the probe 
field, (b) Three-level cascade configuration. To obtain signif- 
icant nonlinear effects for the probe field with Rabi frequency 
LJi, a large detuning Si from the intermediate level |2) is nec- 
essary. 



ular system in the presence of EIT, usually arise by one 
of the two following related mechanisms. One is to vi- 
olate the strict two-photon resonance condition, with 
a frequency mismatch smaller than the width of the 
transparency window [3, Il4l . Alternatively, one can 
add additional energy level(s) in order to induce an ac- 
Stark shift and effectively tune the signal out of reso- 
nance J), 10, 16] . Both mechanisms result in large nonlin- 
earities, accompanied by very weak absorption. Recently, 
the so-called M-scheme, shown in Fig.|21 has been studied 
and proposed as a promising source of giant nonlineari- 
ties that can be utihzed for XPM [ll[i3- The double A 
nature of this M configuration offers the opportunity of 
a simultaneous group velocity reduction for pulses prop- 
agating inside the atomic sample. Group velocity match- 
ing, originally pointed out by Lukin and Imamoglu |l8l |. 
is important to obtain a large XPM. In fact, it has been 
shown by Harris and Hau that if equal group velocity 
reduction is not achieved for both fields, the nonlinear 
phase accumulation will saturate at a certain constant 
value. The consequence is that increasing the length of 
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the sample in which the nonhnear interaction takes place 
is not useful. On the other hand, if group velocities are 
equal, the nonlinear phase accumulation becomes linear 
in the interaction length and it may become very 
large. 

A large cross Kerr phase shift is very useful for 
photonic-based implementations of quantum information 
(QI) processing systems 4, 19]. In fact, a fundamental 
building block for quantum information processing is the 
quantum phase gate (QPG). In a QPG, one qubit gets 
a phase conditional to the other qubit state according to 
the transformation 0, K)i|j)2 — > exp |z)i|j)2 
where {i,j} = 0, 1 denote the logical qubit bases. This 
gate is universal when the conditional phase shift (CPS) 

= 011 + 000-010-001, (1) 

is nonzero, and it is equivalent to a CNOT gate up to 
local unitary transformations when = tt 0, l20j . 

To obtain a CPS of = tt, one looks for a strong in- 
teraction between qubits, ideally accompanied by weak 
decoherence. Photons are a particularly attractive choice 
for qubits due to their robustness against decoherence 
during the processing and transmission of information. 
This feature should ideally permit the transmission of the 
quantum information stored in very weak quantum pulses 
over very long distances with a negligibly small reduction 
of the initial signal. There is however an important dif- 
ficulty in the implementation of an all-optical-QPG: to 
process the information one needs strong photon-photon 
interaction. In fact, to implement QI with photons, a 
nonlinear interaction is needed either to build a two- 
photon gate operation 0, 0, |^ or at the detec- 
tion stage in linear optics quantum computation It 
should also be mentioned that the generation of single- 
photons (which is also necessary in linear optics quantum 
computation) also relies on nonlinear interactions. 

In this paper we perform a semiclassical analysis of 
the interaction of light with atoms in the M configura- 
tion, in which the amplitude of the four fields involved is 
described in terms of the corresponding Rabi frequency. 
The aim is to estimate the effects of noise sources, such 
as dephasing and spontaneous emission, both on the non- 
linear interaction and on group velocity matching. The 
semiclassical regime offers a clear picture of the physical 
aspects involved in EIT-based nonlinear opt ics, and well 
describes a number of recent experiments [Tsl . To 
this end we consider two different configurations of atom- 
field interactions, which we will call the asymmetric (see 
Fig. 121) and the symmetric (see Fig. I1U|I M scheme. The 
paper is thus composed of two main parts. In Sec. ^ 
we describe the physics of the asymmetric M-schemc. 
We start by defining the system and calculating the sus- 
ceptibilities using an approximate treatment employing 
amplitude equations. These analytical calculations are 
then compared with the results of the numerical solution 
of the full system of Bloch equations. Finally, the condi- 
tions for group velocity matching are analyzed. In Sec. lIIII 



the physics of the symmetric M-scheme is described by 
following the same order as in Sec. ^ Conclusions are 
drawn in Sec. IIVI 

II. THE ASYMMETRIC M SCHEME 

A. The System 

The Af-system under consideration has a double adja- 
cent A structure as shown in Fig.|21 where atoms with five 
levels (three ground states |1), |3), |5), and two excited 
states |2), |4)) interact with four electromagnetic fields. 
This configuration can be realized in Zeeman-splitted al- 
kali atoms, such as ®^Rb atoms. The Rabi frequencies 
associated with the lasers driving the atomic transitions 
are defined as 

^ (2) 

where £k is the electric field amplitude, /ly is the rel- 
ative dipole matrix elements induced on the transition 
N) b')- Oil transitions |3) ^ |2) and |5) ^ |4) we 
apply two strong fields, the coupler Q2 and the tuner 
respectively. On the transition |1) <-> |2) a probe field 
is applied (with fii), while on the transition |3) ^ |4) 
a trigger field (with rt^) is applied. In this paper, we 
will analyze the XPM and the group velocity matching 
between the probe and the trigger fields. We call the 
scheme of Fig. 2 the asymmetric M scheme due to the 
asymmetric distribution of the initial atomic population. 
All the atoms are in fact assumed to be initially in state 
|1) so that they directly feel the effect of the probe field 
only, while the effect of the trigger field is only indirect. 
Due to this inherent asymmetry, the dynamics experi- 
enced by probe and trigger fields are always different, 
even when the corresponding parameters (Rabi frequen- 
cies, decay rates, detunings) are equal. The symmetric 
version of this scheme will be analyzed in Sec. IIIII 
The detunings Si (see Fig. ^ are defined as follows 



E2-E1 = TiLOi + hSi (3a) 

E2-E-i = huj2 + hS2 (3b) 

Ei-Es = hL03 + hS3 (3c) 

E4 — i?5 = huj4 + hSi, (3d) 



where (i = 1, . . . , 5) is the energy of level and toi 
is the frequency of the field with Rabi frequency Sl^. 
The Hamiltonian of the system is 

5 

Ha = ^i?.|z)(i|+;i(r!ie-^'^i*|2)(l| + r!2e~^'^=*|2)(3| 

i 

+ r!3e-*'"^*|4)(3| + r!4e"*'"^*|4)(5| + h.c.) , (4) 
where h.c. denotes the hermitian conjugate. Moving to 
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FIG. 2: Asymmetric M scheme. The probe and the trig- 
ger fields, with Rabi frequencies Qi and Q3 respectively, to- 
gether with the stronger pump fields, the coupler and the 
tuner (with Rabi frequencies Q2 and Q4, respectively) drive 
the corresponding transitions. All the atoms are assumed to 
be in state |1) and the detunings are defined in Eqs. (j^J. 



the interaction picture with respect to the following free 
Hamiltonian 



Ho - Ei\l)(l\ + {E2-hSi)\2)(2\ 
+ (£;4-Mi3)|4)(4| + (£;5- 



iE3-hSi2)\3){S\ 



Mi4)|5)(5|, 



(5) 



(6a) 
(6b) 
(6c) 



where 

512 = Si - S2, 

513 = 61-62+ 63, 
6ii = 61 - 62 + 63 - 64, 

we get the following effective Hamiltonian 

H^ff = Mi|2)(2| +Mi2|3)(3| +Mi3|4)(4| +Mi4|5)(5| 

+hni\2){i\ + hn2\2){3\ + hn3\4:){3\ + hn^\4:){5\ 
+nni\i){2\ + fLn*\3){2\ + hni\3){4\ + nf]*|5)(4|. (7) 

B. Amplitude variables approach 

We now study the dynamics driven by Eq. I?)). How- 
ever, we have to include the effects of spontaneous emis- 
sion and dephasing, and we first treat them in a phe- 
nomenological manner by including decay rates Tf^ for 
each atomic level \i) in the equations for the amplitude 
variables (AV) of the atomic wave-function. From an in- 
tuitive point of view, for the excited levels |2) and |4) 
these rates describe the total spontaneous decay rates, 
while for the ground states the associated decay rates de- 
scribe dephasing processes Therefore, the evolution 
equations for the amplitudes bi (t) of the atomic state 



(8) 



become 

bi = 
62 = 

k = 

h = 



AV 



-hi — ifllb2, 



^ 2 



rAV 
^ 3 



rAV 
J- 4 



rAV 
^ 5 



(9a) 

+ i6ijb2-inibi-in2b3, (9b) 
+ i6i2 ] &3 - if^2^2 - i^lbi, (9c) 



i(5i3 j 64 - ifi3&3 - ^^1465, (9d) 
i6u I 65 - i^4b4. (9e) 



The system's initial state is assumed to be the ground 
state Since an efhcient XPM requires a dispersive in- 
teraction, we tailor the dynamics in such a way that this 
initial condition on the populations remains essentially 
unaltered, even when the system reaches the steady- 
state, i.e.. 



1. 



(10) 



To this end we assume that the control field fl2 is stronger 
then the probe field Qi, with the system being approxi- 
mately on Raman resonance for the first and the second 
A subsystems (^i ~ 62 and ^3 ^ 64). Equations © are 
then solved in the steady-state. In order to get a consis- 
tent expression for the nonlinear susceptibilities one has 
to consider higher order contributions to Eq. (|1U|I . which 
is obtained by imposing the normalization of the atomic 
wave-function of Eq. © at second order in |r2i/ri2|. One 
gets the following expression for the steady state ampli- 
tudes 
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where we have defined 



d2 


= 6i- 


-^Vf|2^ 


d3 


= 612 




d4 


= 613 




d5 


= 614 





D„ 



\d2d3 



\n2?\ 



\d4dz 



If^4pl 



d2d^\n3?. 



(11a) 

(lib) 
(11c) 
(lid) 
(lie) 



(12a) 
(12b) 
(12c) 
(12d) 



(13) 
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These results can be used to determine the probe and 
trigger susceptibiUties, which are defined as 



XP 

XT 



A^A*34 ,ss,ss,* 
-O4 O3 



(14a) 
(14b) 



where N is the number of atoms interacting with the elec- 
tromagnetic field, V is the volume occupied by the gas, 
and eo is the vacuum dielectric constant. Doppler broad- 
ening is neglected here. It is well known that first order 
Doppler effect can be cancelled by using co-propagating 
laser fields ,6:j. In particular we emphasize that this is 
valid for cold atomic media in a magneto-optical trap as 
well as for a standard gas cell. 

Inserting Eqs. l(TT|l into Eqs. l(Tl|l and expanding in 
series at the lowest orders in the probe and trigger electric 
fields, f 1 and £3 respectively, one gets 



XP 

XT 



x'P 

(3,cfe) 



(3,cfe) 



(15a) 
(15b) 

where we have introduced the linear susceptibility Xp'', 
the third-order self-Kerr susceptibility Xp'*'^'' ^-iid the 



third-order cross-Kerr susceptibilities Xpt'^'- Eqs. H15|l 
clearly show the asymmetry of the scheme between the 
probe and trigger fields, with the latter possessing a 
nonzero cross-Kerr susceptibility only. This is a conse- 
quence of the asymmetry of the population distribution, 
which essentially remains in the ground state |1) all the 
time. This means that the trigger field drives a virtually 
empty transition, hence the contribution to the suscepti- 
bility comes only from higher order (see for discus- 
sion on the link between the population distribution and 
a linear contribution to susceptibility). It will be shown 
in Sec. IIIII that the symmetric M-scheme brings about 
both a linear and a self-Kerr contribution to the trigger 
susceptibility. 

By using Eqs. (|ll|l and the definitions of Eqs. into 
Eqs. (|14|l . and comparing with Eqs. H15I) at the corre- 
sponding order in the electric fields, one gets the explicit 
dependence of the linear and nonlinear susceptibilities as 
a function of the system parameters, i.e., 



N\^i 



121 



312 



zrr/2 



Ffeo (^1 - (<5i2 - iTi^l2) - l^iaP 

(16) 

for the probe linear susceptibility, and 





\6r. 










(<5i-zrf^/2) (<5i2-«r^^/2)-|r!2P 


2 



Xp 



{3.ck) _ ^|M12P|M34P 

Xp — 



|r!2|2(5i4-ir^V2) 



{3,ck) A^|/^12p|/U34p 
Xt = 



vh'eo [{s, - (<5i2 - *r^^/2) - \n,\^Y [{6,3 - (<5i4 - *r^^/2) - |r!4r 

|f^2p('5i4-ir^^/2) 



(17a) 



Vh so \{S, - lYf'll) (<5i2 - iVfll) - W\ (<5i3 - *rf^/2) (<5i4 - iVf'll) - 



-(17c) 



for the third-order nonlinear susceptibilities. The 
two cross-Kerr susceptibilities are identical whenever 



typical EIT situation we are considering in which 
is large enough. In fact, when ^ 



the quantity (<5i - iVfji) ((5i2 - ir^^/2) - is - ir^^/2) (^12 - ir^^/2)|, one has [T| 

(at least approximately) real. This happens in the 



(3,cfe) 
Xp 



(3,cfc) 
AT 



Jl4 



/2 



\^2? (-513 - ^/2) (<5i4 - ivf'li) ~ |i74r 



(18) 



We shall see in the next subsection that these approxi- of the system. 

mate expressions for the nonlinear susceptibilities fit very The asymmetric M-scheme can be seen as an extension 
well with the numerical solution of the exact dynamics of the four level A^-scheme introduced in Ref. Q , with the 
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addition of the coupling to an additional level |5) pro- 
vided by the tuner field with Rabi frequency f24. In fact, 
it easy to check that upon setting = in Eq. (|17b ). 
one recovers the third-order nonlinear susceptibility of 
the four-level A^-scheme derived in Refs. Hill). As we 
will see below, the role of the tuner field is to enable a fine 
tuning of the group velocities, in order to achieve group 
velocity matching between probe and trigger | 14il6i.il8 ]. 

C. Comparison with the Optical Bloch Equations 

We now study the dynamics of the asymmetric M 
scheme of Fig. |5] by means of the optical Bloch equa- 



tions (OBE), which allow to describe spontaneous emis- 
sion and dephasing rigorously and no more phenomeno- 
logically as in the AV treatment presented in the preced- 
ing subsection. We consider six spontaneous decay chan- 
nels, i.e., the decay of the excited state |2) onto the three 
ground state sublevels |1), |3), and |5) with rates r2i, 
and respectively, and the corresponding decay of 
the excited state |4) onto the three sublevels |3), and 
1 5) with rates r4i, r43 and r45 respectively. Moreover 
we consider dephasing of the each level |z) with dephas- 
ing rate so that the master equation for the atomic 
density operator p is given by 



r ^ 

[Hts% P] + X! X! 'Y {^^^^P^kl - ^ll^klP - P^il^kl^ + X! i'i^kkP^kk - CTkkP - PCTkk) , (19) 

i=2,4 /c=l,3,5 k=l 



where H^^^^ is given by Eq. {Tj) and (iki = The cor- 

responding system of OBE's for the mean values (Jij (t) = 
{<Jij{t)) = Pji{t) is displayed in Appendix 1X1 as Eqs. (|A1|) 
and ljA2p , where we have defined for convenience the total 
decay rates 

r2 = r2i + r23 + r25, (20) 
r4 = r4i + r43 + r45, (21) 

and the composite dephasing rates 

i = 1...5. (22) 

The OBE for the M scheme are quite involved and less 
suited for an approximate analytical treatment with re- 
spect to the AV equations of the preceding subsection. 
In fact, if we consider again the condition |f7i/f22| 1 
and, consistently with Eq. (|10|l . we assume that 

an « 1, (23a) 
a,, « 0, J =2,..., 5, (23b) 

at the steady state, it is possible to see that by insert- 
ing Eqs. (|23|l into Eqs. (|A2I) for the coherences, one gets 
a satisfactory expression for the probe linear suscepti- 
bility only. To be more specific, only the approximate 
linear susceptibility fits well with the numerical solution 
of the OBE, while it turns out to be extremely difficult 
to derive analytical expressions from Eqs. (|A1|) and (|A2p 
for the nonlinear susceptibilities, as simple as those of 
Eqs. IjlTI) . and which reproduce in the same way the ex- 
act numerical solution in the EIT regime we are studying. 
Obviously, one can exactly solve analytically the OBE, 
but the resulting expressions are very cumbersome and 
not physically transparent such as those of Eqs. (|17|l . For 



this reason we will analytically derive from the OBE the 
probe linear susceptibility only, and we will then use the 
OBE only for the numerical determination of the atomic 
steady state. Additionally, deriving this result will en- 
able us to draw a formal analogy between the AV and 
OBE treatments (see Eqs. below). 

The probe susceptibility is defined in terms of the 
atomic coherence cri2 as (see also Eq. (E^)) 

A^M12 N\iXi2? 
XP = p gl2 = „ cri2. (24) 

Using Eqs. (|23|) and performing a series expansion at the 
lowest order in the probe and trigger fields, we arrive 
at an approximate solution for (T12, which, inserted into 
Eq. (|24|l . gives the following expression for the probe lin- 
ear susceptibility 

(1) _ A^|Mi2p 5i2 - ni3/2 

" vTiEQ [<5i2 ~ 1713/2] [<5i- 1 (r2 + 712) /2]-|r!2P' 

. (25) 

By comparing Eq. H25|) with Eq. (|16|) . one can imme- 
diately see that the AV and OBE predictions for the 
probe linear susceptibility coincide provided that the 
phenomenological decay rates Vf^ are appropriately in- 
terpreted, i.e., 

^ r2+7i2, (26a) 
rr ^ 713. (26b) 

This comparison shows therefore that the AV approach 
provides a treatment of the atomic dynamics simpler 
than the OBE's approach, but roughly equivalent, and 
that the intuitive interpretation of its phenomenologi- 
cal decay rates as spontaneous emission total decay 
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rates for the excited states, and as dephasing rates in the 
case of ground state sublevels, is essentiahy correct, es- 
pecially in the typical case in which dephasing rates are 
much smaller than spontaneous emission decay rates (see 
Eqs. (ES). 

We then consider the numerical solution of the OBE 
and we compare it with the analytical treatment based 
on the AV approach presented above. The numerical 
calculations are performed in the range of parameters 
corresponding to KIT, i.e., |rii|, |ri2| ^ l^^sL |^^4| and we 
stay near two-photon resonance for both the probe and 
the trigger field. In Figs. I3lb1 we compare the analytical 
solutions of Eq. IjlGfl and Eqs. H17() with the numerical 
solution of the complete set of Bloch equations given in 
the Appendix A. From these plots it is evident that our 
analytical treatment works satisfactorily well, except for 
a small interval of values of the detuning, corresponding 
to the maximum probe (or trigger) absorption. In such 
a case, the detunings match the Rabi frequencies of the 
two pumps, and the probe (or trigger) field is in reso- 
nance with a single atomic transition. The atoms are 
significantly pumped to the excited levels and the pop- 
ulation assumption of Eq. (|10|l is not fulfilled. In fact, 
the discrepancy between the exact numerical solution of 
the OBE and the AV approach is strictly related to the 
atomic population out of level |1) which, in the case of 
Fig. 13 is about 14% of the total population. 

Figs. l3l6l refer to a situation with small dephasing rates 
(V i,j, 7y = F^^ = F^^ = 10-''r4 ~ few kHz) which are 
typical for not too dense gases. For larger values of the 
dephasing rates (some tens of kHz), we have seen that the 
analytical prediction of the AV approach of the preceding 
subsection starts to depart from the exact solution of the 
OBE. 




FIG. 3: Comparison of the numerical solution (dotted line) of 
the OBE with the analytical prediction of Eq. H16|l (full line) 
for the real part (above) and imaginary part (below) of the lin- 
ear probe susceptibility versus the normalized probe detuning 
Si/Fi. The parameter used are the following: = r2 = 36 
MHz, ri^ = r4 = 38 MHz, (Jz = (S3 = ^4 = 0, V 



02 = 03 = 04 = u, V 1,J 

= ri^ = r^^ = io~^r4, ni = o.08r4, ^2 = 2r4, 



!^3 = 0.04r4 



4 

^^4 = r4. 



so that, using Eqs. p5|l . the nonlinear cross-phase shift 
for the two fields of interest is given by 



6p = ^ / dzRe[x'/']\eT{z,t)\\ 



(30a) 



^t=yJ rf^Re[x?;''] |£p(^,t)|' , (30b) 



D. Group velocity matching 

The propagation equation for the slowly varying elec- 
tric field amplitudes £i{z, t), i = P,T, defined as 

£i{z^ t) = £i{z^ t) exp {ikiZ — icuit} + c.c. i — P^T, 

is given by 

d 1 d \ k 
oz Vg ot J 2 

(27) 

where Vg is the group velocity, generally defined as Vg = 
c/(l -|- rig), with c the speed of light in vacuum and 



(28) 



where I is the length of the atomic medium. These non- 
linear cross-phase shifts are of fundamental importance 
also for quantum information processing applications. In 
fact, the CPS of Eq. is determined only by these cross- 
Kerr contributions to the total phase shift, because the 
linear and self-Kerr contributions cancel out, as shown in 
Refs. 14, 15J. 

For Gaussian probe and trigger pulses of time dura- 
tions Tp and tt, and with peak Rabi frequencies ^Vp°'^ 
and VL^°'^ respectively, the nonlinear cross-phase shifts 
can be written as (see also Refs. [T3. iTsj l') 



4c 



])rp 

IM12I 



M34r 
2|r-vPeafe|2 



4c 



Cp 
erfjCx' 
Ct 



Mxp J 



(31a) 
(31b) 



the group index, oji being the frequency of field i. The 
solution of Eq. (|77ll is 

e,(z,<)=£,(0,i-^)exp|z.|^ dz'x.(z', t)| , (29) 



where Cp = {1 — Vg /vJ)^/2l/VgTT and Ct is obtained 
from Cp upon interchanging the indices P ^ T . Large 
nonlinear cross-phase shifts take place for appreciably 
large values of the two cross-Kerr susceptibilities real 
parts, and especially when probe and trigger velocities 
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FIG. 4: Comparison of the numerical solution (dotted line) 
of the OBE with the analytical prediction of Eq. llTb 'l (full 
line) for the real part (above) and imaginary part (below) 
of the probe cross-Kerr susceptibility versus the normalized 
probe detuning S1/T4. To reduce as much as possible the 
influence of the self-Kerr susceptibility we have considered a 
probe Rabi frequency f2i much smaller than that of the trigger 
field. Parameters are: T^^ = = 36 MHz, ^ = r4 = 38 

MHz, 52 = (53 = 54 = 0, V i,j 7,, = r^*^ = r^^ = io-'*r4, 
0.1 = o.oo4r4, ^2 = 2r4, = o.o4r4, ^4 = r4. 



FIG. 5: Comparison of the numerical solution (dotted line) 
of the OBE with the analytical prediction of Eq. 1117b ') (full 
line) for the real part (above) and imaginary part (below) of 
the probe self-Kerr susceptibility versus the normalized probe 
detuning 5i/r4. To reduce as much as possible the influence 
of the cross-Kerr susceptibility we have considered a trigger 
Rabi frequency fia much smaller than that of the probe fleld. 
Parameters are: T^^ = Ta = 36 MHz, Tf^ = r4 = 38 



MHz, 52 = 53 = 54 = 0, V i,j 7,, = r 

fii = o.5r4, ^2 = 2r4, n-i = o.oo5r4 



3 — i 5 

= r4. 



= io-*r4 



where ll^ 



become equal, i.e., when Cp,t — > 0, in which case the 
erf[C]/C reaches the maximum value In this limit 

the cross-phase phase shifts linearly increase with the 
length of the atomic medium /. This explains why achiev- 
ing group velocity matching, Vg = wj, is of fundamental 
importance. Moreover group velocities become small for 
large group indices and this condition can be achieved 
within the EIT transparency window, where Re[x] van- 
ishes, and the group velocity is strongly reduced due to 
a large dispersion gradient dRe[x]/dLu. 

Let us see how small and equal probe and trigger group 
velocities can be obtained. We consider the approximate 
analytical expressions for the susceptibilities of Eqs. H15(l - 
(|17|l derived above within the AV approach, and which 
we have seen to work very well in the EIT regime. As- 
suming to stay at the center of the transparency window 
for the probe {S12 — 0) where the dispersion gradient 
is maximum, and neglecting dephasing rates and 
F^^, which are typically much smaller than all the other 
parameters, one gets 
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V 2?ieo|f^2p 



(1 



(32a) 
(32b) 



f3 



tuation we 
so that, using Eqs. 



(33) 

In the EIT situation we are considering it is 71^, ^ 1, 



c 
c 



2heoc\n2\^ 



{N/V)\fii2\^uJi{l + \ns\^P)-' 
2heoc\n2\^ 



nT - (iV/t/)|M34Pc.3|OiP/3- 



(34a) 



(34b) 



As expected, the asymmetric M-scheme does not yield 
equal slow down of both trigger and probe pulse automat- 
icall y as , for example, the scheme of Petrosyan and Kur- 
izki |l6j does. In fact, the two expressions of the group 
velocities are generally different. Nonetheless, Eqs. H34|l 
show that group velocity matching is always achievable 
by properly adjusting the parameter (3, which means ad- 
justing the tuner intensity |r24p and the composite de- 
tuning (5i4. This shows that the present asymmetric M- 
scheme can be seen as a modified version of the iV-scheme 
of Ref. ||5j, in which the tuner pump field is added just 
in order to "tune" the group velocity of the trigger pulse 
so to make it equal to that of the probe. The possibility 
to achieve group velocity matching is shown in Fig. 
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FIG. 6: Comparison of the numerical solution (dotted line) 
of the OBE with the analytical prediction of Eq. il7t ) (full 
line) for the real part (above) and imaginary part (below) 
of the trigger cross-Kerr susceptibility as a function of the 
normalized trigger's detuning Ss/Ti. Parameters are similar 
to those of Fig. H = Fa = 36 MHz, Ff ^ = F4 = 38 

MHz, 5i = (52 = 54 = 0, V i,j = F^^ = F^^ = 10-*F4, 
Qi = O.O8F4, ^2 = 2F4, ^3 = 0.04F4, ^4 = F4. 



where both the numerical result derived form the OBE 
and the approximate analytical expressions of Eqs. H34|) 
are plotted versus the trigger detuning ^3. Two different 
values of ^3 exist for which = ~ 1000 m/s (see 
Fig. d) . Parameter values here correspond to typical val- 
ues for a cell of ^^Rb atoms, i.e., T^^ = Fa ~ 36 MHz, 
rf ^ = r4 ~ 38 MHz, N/V ~ 3 X lO^^ cm^^^ 61=62 = 0, 

(54 ~ (53 - 2or4, fii = o.08r4, = 2r4, = omt^, 

= r4, V i,j, 7y = = = 10-4r4. More- 

over Fig. [3 clearly shows that the simple expressions of 
Eqs. (|34() well reproduce the exact numerical solution of 
the OBE. 



E. Pulse propagation 

In previous Section, we have addressed the problem 
of group velocity matching between probe and trigger 
fields in the asymmetric M-scheme. It should be empha- 
sized that the analysis and the results presented there 
are strictly valid for the continuous-wave (cw) fields. 
We would now address that same problem but with the 
pulsed probe and trigger fields in mind. At the first 
look, Eqs. (Ell appear to suggest that the group velocity 
matching would not be possible in the pulsed regime. As 
the group velocity of the trigger pulse is inversely pro- 
portional to the square of the probe pulse, trigger suffers 
anomalous dispersion, i.e. in the presence of a pulsed 
probe, the trigger pulse will get distorted, splitting into 
several components, each having a different group veloc- 




FIG. 7: Group velocity of the probe and trigger pulses versus 
the normalized trigger detuning 53/F4. Full lines denote the 
analytical predictions of Eqs. (1341 (the thick line refers to 
the probe and the thin line to the trigger). Circles and dots 
refer to the numerical solution of the OBE for the probe and 
trigger group velocity, respectively. This figure shows how it is 
possible to obtain group velocity matching in the asymmetric 
M-scheme: two different values of 63 exist for which = 
vj ~ 1000 m/s. The parameters are those of the Di and D2 
line in the *^Rb spectrum: F^^ = Fa ~ 36 MHz, Ff ^ = 
F4 ~ 38 MHz, 5i = 62 = 0, 54 ^ (53 ~ 2OF4, ^1 = O.O8F4, 
Q2 = 2F4, fia = 0.04F4, = F4, V i,j, -/ij = Ff^ = Ff^ = 
W"^r4, N/V = 3.0 ■ 10^^ cm-^ 



ity. 

It will be shown in this Section that the above con- 
clusion is an artifact of approximations made to obtain 
a closed and compact expression for group velocities. In 
particular, the pulse propagation in this approximation 
is described by Eqs. H27I) . with group velocities given 
by Eqs. 1)34(1 and nonlinear susceptibilities Xi being those 
of Eqs. H17|l . This is equivalent to the adiabatic elimi- 
nation of the atomic degrees of freedom. Such adiabatic 
elimination, strictly speaking, is not valid in the param- 
eter regime explored in this paper: strong nonlinear in- 
teraction between probe and trigger pulses suggests that 
the contribution of the atomic medium is far from being 
adiabatic. Also, it should be noted that the dephasing 
processes have been neglected in the derivation of Vg. For 
the adiabatic case, the above conclusion is correct: the 
trigger pulse suffers anomalous dispersion and its group 
velocity becomes singular towards the edges of a probe 
pulse. However, pulse propagation through the asym- 
metric M-system do not follow such a simple approximate 
evolution. Full propagation problem must then be solved 
which includes adding the time-dependent equations for 
the pulses 
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probe pulse trigger pulse 




time [a. u.] ^ -50 z [a. u.] time [a. u.] -50 ^ [a. u.] 



FIG. 8: (Color online) Propagation of probe and trigger pulses through the asymmetric M medium. Pulses are taken to be 
Gaussian at time t = and are sufficiently long (r^ > i = P,T. Units are arbitrary, with c = 1. 



probe pulse 




time [a. u.] ° '^^ z [a- u.] 



FIG. 9: (Color online) Propagation of probe and trigger puL 
Gaussian at time t = and short (r^ < l/AtjJ,.), i — P,T. Un 



trigger pulse 



0.04r 




time [a. u.] " z [a. u.] 



through the asymmetric M medium. Pulses are taken to be 
are arbitrary, with c = 1. 



to the OBEs [Eqs. HAlllA2p ]. and numerically solving the 
resulting system of equations. In the above equation, it is 
understood that ap — a 12 and gt — C34 in the notation 
of Eqs. (IX2I1 . 

Results are shown in Figs. |S1 and El for the same set 
of parameters that yields group velocities matching in 
Fig. □ Vertical axes have been scaled appropriately to 
obtain Rabi frequencies fi^. Two operating regimes could 
be identified, long-pulse regime and short-pulse regime, 
where 'long' and 'short' denotes the pulses' length in 
time. This length is compared to the inverse width of 
the transparency window. Long pulses fit well into the 
transparency window, while short pulses do not. Fig. |S1 
shows the results of our simulation for the initially identi- 
cal long Gaussian pulses. It is clear that the pulses prop- 
agate undistorted with the equal group velocities. Tiny 
amplitude decay is present due to the small imaginary 
part of the nonlinear susceptibility. 

Short pulses (Fig. |5J| however, show distortion. Probe 



pulse distortion comes from the absorption, as the pulse 
spreads outside of the transparency window. Trigger 
pulse shows the same absorption effect, but moreover it 
also splits into several components which then continue 
to propagate with a different group velocities each. Note 
that the singularity present in the adiabatic approach is 
not present here. This is due to the fact that the de- 
phasing, neglected in the adiabatic treatment, effectively 
regularizes the equations, removing the singularity. 



It is also noted that in the long-pulse regime, both of 
the pulses propagate virtually undistorted, with a group 
velocity uniform across each of the pulses. Our simula- 
tions suggest that the approximate Eqs. are valid, 
as long as the Rabi frequencies there are considered to 
be taken at the peak of the pulse, i.e. 17^ — > Vt^""^ ■ 
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III. THE SYMMETRIC M-SCHEME 

In this section we analyze the symmetric M~scheme, 
schematically shown in Fig. ^| The initial conditions 
and the configuration of the fields are slightly different 
from those of the asymmetric case of Sec. ^ The same 
five levels could be used, but all the atoms are now ini- 
tially prepared in level |3) (see Fig. I10|l . Moreover, the 
role of the probe and of the coupler fields are exchanged, 
i.e., now the probe field (still with Rabi frequency Qi 
and central frequency oji) couples levels |2) and |3), while 
the coupler (still with Rabi frequency and central fre- 
quency UJ2) induces transitions between levels |1) and |2). 
The role of trigger and tuner fields remains unchanged. 
In such a way, the scheme becomes symmetric for probe 
and trigger, and the two fields experience exactly the 
same dynamics whenever the corresponding parameters 
are made equal, i.e., when the Rabi frequencies are cor- 
respondingly equal (Oi = fls, = fli), as well as the 
detunings, {Si — 63, S2 = Si), which are now defined sim- 
ilarly to those of the asymmetric M scheme (see Eqs. ©) 



except for probe-coupler exchange, i.e., 

E2-E1 = hLU2 + hS2, (36a) 

E2- E3 = huji + hSi, (36b) 

Ei-E3 = hLU3 + hS3, (36c) 

E4-E5 = huj4 + nS4. (36d) 



In this way, the scheme can be seen again as formed 
by two adjacent A, one for the probe and one for the 
trigger, now however symmetrically placed with respect 
to state 1 3). As we have done for the asymmetric M 
scheme, we assume to stay close to the two-photon res- 
onance conditions. Si ~ S2 and ^3 ~ ^4, and more- 
over that |fii| <C 1^2], and jfisl <C \^4\, so that both 
probe and trigger will experience FIT. As we have seen 
above, a large XPM is obtained when the group velocities 
are equal [la. El El El El El , and the advantage of 
the present symmetric M-scheme is that group velocity 
matching is automatically achieved once that the scheme 
is exactly symmetric between probe and trigger. 
The Hamiltonian of the system is 

5 

Hs = ^i?,|z)(z| + ?i(r!ie-^"^*|2)(3|-t-r!2e-'-^*|2)(l| 

i 

+ f73e-*"'^*|4)(3| + f74e-*'^-'*|4)(5| + h.c.) . (37) 

Moving to the interaction picture with respect to the 
following free Hamiltonian 

i/^ - E3\S){3\ + iE2-hSi)\2){2\ + iEi-hSi2)\l){l\ 
+ iEi - M3)|4)(4| + (^5 - M34)|5)(5|, (38) 

where 

S12 = Si - (52, (39a) 
S34 = S3-S4, (39b) 



4) 




FIG. 10: Symmetric M scheme. The probe and the trig- 
ger fields, with Rabi frequencies Oi and Q3 respectively, to- 
gether with the stronger pump fields, the coupler and the 
tuner (with Rabi frequencies ^2 and SI4, respectively) drive 
the corresponding transitions. AH the atoms are assumed to 
be in state |3) and the detunings are defined in Eqs. I36II . 



we get the following effective Hamiltonian 

H^ff = Mi|2)(2| + Mi2|l)(l| + M3|4)(4| + M34|5)(5| 

+nni\2){3\ + hn2\2){i\ + hn3\4){3\ + ?if]4|4)(5| 

+Kll\3){2\ + hn*2\l){2\ + nn*\3){4\ + hni\5){A\. (40) 

A. Amplitude Variables Approach 

We first study the system dynamics by means of the 
AV approach, in which the state of the atom is described 
by the wave- function of Fq. whose time evolution 
is determined by the Hamiltonian of Eq. (|4U|I . supple- 
mented with phenomenological decay rates F^^ for each 
atomic level The corresponding evolution equations 
for the amplitudes bi{t) are 



bi{t) - 


—idibiit) 


-in*b2{t), 


(41a) 


kit) = 


-id2&2(i) 


- in2biit) - iflibsit), 


(41b) 


kit) - 


-1^363 - 


in{b2it) - in^biit), 


(41c) 


kit) = 


—id^biit) 


- iQ.3b3it) - lilib^it), 


(41d) 


kit) = 


-4^565 (i) 


- in^biit), 


(41e) 



where, similarly to what we have done for the asymmetric 
case, we have defined 

di = Si2-irt^/2, (42a) 

d2 = Si-iYfl2, (42b) 

d3 = -iVf 12, (42c) 

di = S3-iVfl2, (42d) 

ds = S3i-iYf 12. (42e) 

Since we choose again |f2i/f22| ^ 1 and |f23/r]4| ^ 1, 
it is reasonable to assume that the atomic population 
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remains in the initial state |3) to a good approximation Eqs. (|14ll l 



6- ^ 1. (43) 

The set of equations (|41|l is then solved in the steady- 
state. In order to get a consistent expression for the 
nonlinear susceptibilities one has to consider higher order 
contributions to Eq. H43|l , which is obtained by imposing 
the normalization of the atomic wave-function of Eq. (jH)) 
at second order in |i7i/fi2| and |il3/r24|. One gets the 
following expression for the steady state amplitudes 



Mil 



in. 



2\did2 



2 \d4dQ 



(44a) 
(44b) 
(44c) 
(44d) 
(44e) 



These results can be used to determine the probe and 
trigger susceptibilities, which are now defined as (see 



7 SS _ 


ilidi 


02 — 


did2 - 




n^ds 




d<sdi - 






&r = 


did2-|f^2p 




i^3^4 


&r = 


4^4 - |r!4p 



Xp 

XT 



VeoSi 



^2 "3 



iV|M32p 

N\^^34\' 



04 03 . 



(45a) 
(45b) 



Inserting Eqs. H44|l into Eqs. I|45(l and expanding in se- 
ries at the lowest orders in the probe and trigger electric 
fields, £1 and £3 respectively, one gets 



Xp 

XT 



xf 

x't' 



Xp 

(3 

Xt 



(3,sk) 



(3,sfc)|^ |2 



3I 



Xp \£3 

{3,ck) I r, 

Xt 



(46a) 
(46b) 



where we have again distinguished the third-order self- 
Kerr susceptibilities Xp'f^ from the third-order cross- 

Kerr susceptibilities Xpt^'^- Using Eqs. H44|l and the def- 
initions of Eqs. (|42() . we get the following expressions for 
the linear susceptibilities. 



iV|M32| 



Ji2-zr^^/2 



vnso (5i - ^T^y|2) (5i2 - *r^i^/2)^^ 

7V|M34p <534 -*r^^/2 



(4? 



vhea (^3 - irf ^/2) (^34 - *r^^/2) - [^4 

and the following ones for the nonlinear susceptibilities 



Xp 



(3.sk) 

Xt 



xy 



(3,cfc) 

Xt 



iV|/i32|^ 






- (<5i2 - ^/2) 


|5i2-*rf^/2f + |r!2|2' 








-^Tr|2) 


(<5i2 - »rf ^/2) - |f^2 


2]|(5i-zr^V2) (<5l2- 


zrf^/2) - |r!2P|' 


^|M34|^ 






-(<534 -^rr/2) 


|534-*r^^/2|' + |r!4|2' 








- ^/2) 


(534 - «r^^/2) - |r!4 


2] |(53-zrf^/2) (<534- 


zr^^/2) - |r!4|2|' 



iV|M32p|M34p 






(^12 - ir^i^/2) 


|(534-^^^^/2|' + |^^4|2 










- iTt''l2) - \n2 


2] |(53-ir^V-/2) (^3^_ 


-zr^^/2) 


-|f^4P|' 


A^|M32P|M34P 






(534-*rr/2) 


|(5i2-zrf^/2|' + |f]2p' 








[(^3 


- zrf ^/2) (^34 




2] |(^i-zr^V2) (<5i2- 


^iT^^/2) 





(48a) 
(48b) 
■,(48c) 
r(48d) 



First of all we note that the expressions of the probe 
and trigger susceptibilities above are completely sym- 
metric. This means that probe and trigger experience 
the same linear and Kerr susceptibilities, as soon as the 
corresponding parameters correspond, i.e., /i32 = /i34, 

= n^, 02 = f^4, h = S3, 62 = S4, r^^ = r^^, 

Tf^ = rf^. Moreover, the probe linear susceptibility 
of Eq. l|T7k ^ and the self-Kerr susceptibility of Eq. |@Ht) 



coincide with the corresponding ones of the asymmet- 
ric case, Eq. (fTC|l and Eq. ifTTk ) respectively, because 
the phenomenological decay rate here plays just 

the same role of the phenomenological decay rate F^^ 
of the asymmetric scheme. This is not surprising, since 
the probe response in the absence of the trigger field is 
the same in the two M scheme studied here. Finally the 
cross-Kerr susceptibilities of the two schemes are gener- 
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ally different, both for the probe and the trigger, even 
though they possess a similar structure. The main rel- 
evant difference between the two cross-Kerr susceptibili- 
ties is in the dependence of their real parts upon the de- 
tunings. In fact, in the asymmetric case both real parts 
are proportional to the composite detuning Su = 612+634^ 
(see Eqs. © and H39|l ). so that one has a nonzero XPM 
as soon as one of the two A subsystem is shifted from 
the two-photon resonance condition. In the symmet- 
ric case instead, Rejxp''^'^''} is proportional to S12 and 

Iie{xx''^'^^} proportional to 634, and the two-photon 
resonance condition has to be violated by both A subsys- 
tems if each field has to experience a nonzero XPM. 



B. Comparison with the Optical Bloch Equations 

We now study the dynamics of the symmetric M 
scheme of Fig. El by means of the OBE, which allow to 
describe spontaneous emission and dephasing more rigor- 
ously. Due to the similarity of the symmetric and asym- 
metric M schemes, we consider the same spontaneous 
emission and dephasing processes described in subsec- 
tion |n3 As a consequence, the master equation for the 
atomic density operator p is again given by Eq. (|19l) , with 
the only difference that the Hamiltonian ^^^j is replaced 
by the corresponding Hamiltonian H^j:j: of the symmet- 
ric scheme, given by Eq. (|40|l . The corresponding system 
of OBE's for the mean values ct^ (i) = (o'ij(t)) = pji{t) is 
displayed in Appendix IbI as Eqs. l|BHl and (|B2(I . where 
we have used the definitions of Eqs. (|2U|) - (|22|I . 

Also in this symmetric case, the OBE are less suited 
for an approximate analytical treatment with respect to 
the AV equations of the preceding subsection. In fact, if 
we consider the conditions |ili/f72| ^ 1 and |il3/fi4| ^ 1 
and, consistently with Eq. (|43|l . we assume that 



fas 



1, 

0, J -1,2,4,5, 



(49a) 
(49b) 



at the steady state, it is possible to see that by insert- 
ing Eqs. (|23() into Eqs. (|B2|) for the coherences, one gets 
a satisfactory expression for the probe linear suscepti- 
bility only. To be more specific, only the approximate 
linear susceptibility fits well with the numerical solution 
of the OBE. It is not easy to derive analytical expressions 
from Eqs. HB1|I and HB2() for the nonlinear susceptibilities 
which would be as simple as those of Eqs. (|48|l and which 
would reproduce in the same way the exact numerical 
solution of the OBE within the EIT regime. Again, one 
could exactly solve analytically the OBE, but the result- 
ing expressions are very cumbersome and not physically 
transparent as those of Eqs. For this reason we will 

analytically derive from the OBE the probe linear sus- 
ceptibility only, and we will then use the OBE only for 
the numerical determination of the atomic steady state. 
In addition, deriving this result will enable us to draw 



a formal analogy between the AV and OBE treatments 
(see Eqs. below). 

The probe and trigger susceptibilities are now defined 

as 



XP 



XT 



Np 



32 



NH34 

Veo£3 



C32 



■0'34 



N\P32? 



0'32, 



-0'34- 



(50a) 
(50b) 



Using Eqs. (|49|l and performing a series expansion at the 
lowest order in the probe and trigger fields, we arrive at 
the following expressions for the probe and trigger linear 
susceptibilities 



(1) ^|M32 

Xp = 



Jl2 



1713/2 



x'^ 



Vheo [612 - 1713/2] [Si 
N\p3a\^ S 



34 



i{T2+7i2)/2]~\&}^^ 



Vheo [S34 ~ ^753 / 2] [S3 -i{T4 + 754) /2] - | 



The probe linear susceptibility of Eq. H51a|l coincides with 
that of the asymmetric M scheme of Eq. (|25ll , as it must 
be, since the linear properties of the probe in the two M 
schemes are identical. Moreover, as noted before, due to 
the symmetry of the scheme, this probe linear susceptibil- 
ity coincides with that of the trigger of Eq. (|51b ) when 
the corresponding parameters coincide, i.e., ^32 — ^34, 
Si = S3, S2 = Si, ^2 = fit, T2 = T4, 754 = 712 and 
753 — 7i3- By comparing Eqs. (|FTTl with Eqs. (|17|l . one 
can also see that the AV and OBE predictions for the 
linear susceptibilities again coincide provided that the 
phenomenological decay rates F^^ are appropriately in- 
terpreted, i.e.. 



T^AV 

^ 2 

rAV 
^ 1 
rAV 
^ 4 
rAV 
^ 5 



^ F2 + 712, 

^ 713, 

^ F4 -t- 754, 

^ 753- 



(52a) 
(52b) 
(52c) 
(52d) 



This shows again that the intuitive interpretation of the 
phenomenological decay rates F^^ as spontaneous emis- 
sion total decay rates for the excited states, and as de- 
phasing rates in the case of ground state sublevels, is 
essentially correct. 

We then consider the numerical solution of the OBE 
and we compare it with the analytical treatment based on 
the AV approach. The numerical calculations are again 
performed in the limits discussed above, i.e., |0i|, [U2[ <C 
[^3[, 1^4 1 and we stay near the Raman resonance for both 
the probe and the trigger. In Figs. llllTl^ we compare the 
analytical solutions of Eqs. (|47|l and Eqs. (|48|) with the 
numerical solution of the complete set of Bloch equations 
given in the Appendix B. Fig. 1111 shows the linear sus- 
ceptibilities and refers to a perfectly symmetric situation 
between probe and trigger, i.e., Tf^ = F2 = Tf^ = 
F4 = F = 27r X 6 MHz, ^1=^3 = 0.08F, = f^4 = F, 
S2 = 54 = 0,V i,j 7,j = Ff^ = F^^ = 10~4f. As 
a consequence, the probe and trigger linear susceptibili- 
ties as a function of the respective detunings Si and S3 
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are two indistinguishable curves. In such a case, group 
velocity matching is automatically guaranteed whenever 
M32 — Fig- E] shows that Eqs. l|T7jl work very well, 
except when the detunings correspond to the maximum 
probe (or trigger) absorption. In such a case, the detun- 
ings match the Rabi frequencies of the two pumps, and 
the probe (or trigger) field is in resonance with a single 
atomic transition. The atoms are significantly pumped 
to the excited levels and the population assumption of 
Eq. 143|l is no more fulfilled. 

Fig. 1121 shows the cross-Kerr susceptibilities again in 
a perfectly symmetric situation between probe and trig- 
ger so that their plots as a function of the respective 
detunings &\ and (^3 exactly coincide. However, in order 
to reduce as much as possible the infiuence due to the 
simultaneous presence of the self-Kerr susceptibility, we 
have considered a probe Rabi frequency 17 1 much smaller 

than that of the trigger in the x'p''^^^ plot and viceversa 
a trigger Rabi frequency fl^ much smaller than that of 
the probe in the x^''^'^' plot. To be more precise. Fig. 1121 
refers to F^^ = F2 = t"^^ 
fi2 = = F, 62 = S4 = 0, \f i,j 7y 
IQ-^'F, and fii = 0.002F, f^g ^ 0.08F in the case of the 
Xp'"''^ plot, and to = 0.002F, fli = 0.08F in the 

case of the Xt '^'^'^ plot. We can see from Fig. E| that 
the AV approach gives a satisfactory description also for 
the cross-Kerr susceptibility, except when probe or trig- 
ger absorption is maximum, as it happens for the linear 
case. 

Finally Fig. shows the self-Kerr susceptibilities 
again in a perfectly symmetric situation between probe 
and trigger. As a consequence their plots versus the re- 
spective detunings 5i and ^3 exactly coincide. Here, in 
order to reduce the influence due to the simultaneous 
presence of the cross-Kerr susceptibility, we have consid- 
ered a trigger Rabi frequency Q3 much smaller than that 

of the probe in the Xp'^''^ plot and viceversa a probe 
Rabi frequency f2i much smaller than that of the trigger 
in the Xt'*'^^ plot. The parameters are the same as in 
Fig. El except that here we have chosen ^2 — = 2F, 



Ff ^ = F4 = F = 27r X 6 MHz, 



rAV 
J- 1 



^ 5 



(3,sfc) 



ni = OAT, Qs = 0.004F in the case of the Xp 

and ^3 = 0.4F, f^i = 0.004F in the case of the Xt'""' 
plot. The agreement between the AV prediction and the 
numerical solution of the OBE is satisfactory. 



C. Group velocity matching 




FIG. 11: Comparison of the numerical solution (dotted line) 
of the OBE with the analytical prediction of Eqs. 11471 (full 
line) for the real part (above) and imaginary part (below) of 
both probe and trigger linear susceptibilities versus their re- 
spective normalized probe detunings Si/T and S3/T. Probe 
and trigger susceptibilities exactly overlap because we con- 
sider the perfectly symmetric situation = r2 = = 
r4 = r = 27r X 6 MHz, ^1=^3= O.OSF, = ^4 = T, 

52 = 54 = 0, V i,j j,j = r^^' = r^^ = 10-^r, /i32 = ^^34, 

which guarantees perfect group velocity matching. 



vantage of the symmetric M scheme over the asymmetric 
one. As we have seen above, the group velocity of a pulse 
is given by 



1 



P T 



(53) 



where the group index rig is given by Eq. H28|l . The con- 
tribution of the nonlinear susceptibilities to Vg is negligi- 
ble with respect to that of the linear one, which is nonzero 
for both probe and trigger in this case (see Eqs. 14t)|) '). 
Therefore, approximating x with the linear contribution 
X*-^-* and inserting Eqs. H47|) into the definition (|28|l . one 
plot, gets the following expressions for the two group indices 



^P _ A^|M32P 

^» 21/neo 

T ^ ^l/^34|' 



Re' 



Rc' 



did2-\n2\^ (dirf2- |f^2pj 



dF,d, 



+ 



{dl + \n.^^ 



5"4 



(^5^4 - |f74p 



§^b,) 



As we have seen in subsection IIIDI the condition of 
group velocity matching is of fundamental importance 
for achieving a large cross-phase modulation between 
probe and trigger fields. It is evident from the inher- 
ent symmetry of the present scheme that the condition 
of equal probe and trigger group velocities is automat- 
ically achieved when the corresponding parameters are 
equal i.e., /i32 = ^J,34, Si = S3, S2 = S4, il2 = ^^4, F2 = F4, 
754 = 712 and 753 = 7i3, uji — W3- This is the main ad- 



where we have used the definitions of Eqs. (|42|l for dj, 
j = 1,2,4,5. The symmetry between probe and trig- 
ger discussed above is evident also in these expressions. 
Eqs. H53I) and 1)54(1 are now compared with the corre- 
sponding ones obtained from the integration of the full 
set of Bloch equations of Appendix ^ The comparison 
is shown in Fig. 1141 which refers to the completely sym- 
metric situation between probe and trigger defined above 
and therefore shows exact group velocity matching for all 
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FIG. 12: Comparison of the numerical solution (dotted line) 
of the OBE with the analytical prediction of Eqs. MSt -d) (full 
line) for the real part (above) and imaginary part (below) 
of both probe and trigger cross-Kerr susceptibilities versus 
their respective normalized probe detunings Si/T and Ss/T. 
Probe and trigger susceptibilities exactly overlap because we 
consider a perfectly symmetric situation: = F2 = F^^ — 
F4 = F = 27r X 6 MHz, ^2=^4= F, ^2 = ^4 = 0, V 
ij 7ij — F^*^ — F^^ — 10~*F; moreover we have chosen 
ill = 0.002F, 0,3 = 0.08F in the case of the Xp''"*' plot, and 
viceversa ils = 0.002F, fli = 0.08F in the case of the x^t '^^'^ 
plot, in order to reduce as much as possible the influence of 
the self-Kerr susceptibilities. 



FIG. 13: Comparison of the numerical solution (dotted line) 
of the OBE with the analytical prediction of Eqs. 1148b . b) (full 
line) for the real part (above) and imaginary part (below) of 
both probe and trigger self-Kerr susceptibilities versus their 
respective normalized probe detunings 5\/T and 83 /T. Probe 
and trigger susceptibilities exactly overlap because we con- 
sider a perfectly symmetric situation: F^^ = F2 = F^^ = 
F4 = F = 27r X 6 MHz, 0.2 = 0.4,^ 2F, (S2 = 54 = 0, V 
iJ 7ij = F^*^ = F^^ = 10~^F; moreover we have chosen 
9.1 = 0.4F, ^3 = 0.004F in the case of the Xp'"''"' plot, and 
viceversa fis = 0.4F, f2i = 0.004F in the case of the Xt'**"' 
plot, in order to reduce as much as possible the influence of 
the cross-Kerr susceptibilities. 



values of the detunings Si =63. Fig. 1141 shows an excel- 
lent agreement between analytical and numerical results. 
The only points in which the two curves do not coincide 
exactly are when the detunings match the Rabi frequen- 
cies of the two pumping field. In fact in this conditions 
the fields are in resonance with a single atomic transition 
and the atoms are pumped to the excited levels. The 
other points that determine the disagreement are in the 
vicinity of the peaks. In fact in these regions the deriva- 
tives are small, because of the change in slope of the real 
part of the susceptibilities. Hence, the group index of 
Eqs. (|54|) is small, and the group velocity jumps near c. 

In some cases, the perfectly symmetric conditions guar- 
anteeing group velocity matching, i.e., /i32 = Si = 
S3, S2 = S4, fl2 = ri4, r2 = r4, 754 = 712 and 753 = 713, 
are difficult to realize in practice. In fact, sometimes 
it may be convenient to use transitions with different 
Clebsch-Gordan coefficients, yielding therefore a signifi- 
cant discrepancy between (132 and (134. The other sym- 
metry conditions above are less problematic because de- 
tunings and Rabi frequencies can always be made equal 
by the experimenter, and moreover decay and dephasing 
rates, even though not perfectly equal, are often compa- 
rable to each other. Just to give an example, one could 
implement the symmetric M scheme of Fig. 1101 by using 
the Di and D2 hue of the ^^Rb spectrum. The Zeeman 



sublevels |5Pi/2F = 1,to = 0) and I5P3/2F = l,m = 0) 
could be chosen as levels |2) and |4), respectively, while 
the Zeeman sublevels \5Si/2F — l,m = —1), \5S1f2F = 
2, m = 1) and \5Pi/2F — 1, m = 1) could chosen as lev- 
els |1), |3) and |5), respectively (see also Ref. for a 
similar choice). For these levels the atomic transitions 
related to the probe and trigger have dipole moment ma- 
trix elements 1132 , /J.34 differing by a factor vTO, violating 
therefore the symmetry condition. It is evident however 
that this slight asymmetry can be compensated (so that 
group velocity matching can be still achieved in a re- 
stricted but still useful range of detunings) by properly 
adjusting the Rabi frequencies of the tuner field and 
of the coupling field f22, which will be no more equal. In 
fact, by imposing group velocity matching at the center 
of the transparency window, i.e., for 612 = (^34 = 0, we 
derive the condition 



fl2 — an4, 
where the correction factor a is given by 



lM32r ^1 

|M34p W3' 



(55) 



(56) 



As shown in Fig. ^1 if the adjustment condition of 
Eqs. (|55|l and (|55|l is taken into account, one still gets 
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FIG. 14: Group velocity of the probe and trigger pulses ver- 
sus the normalized detunings Si/T — S3/T. Lines denote 
the analytical predictions of Eqs. and l|54^ (the full line 
refers to the probe and the dashed line to the trigger) . Circles 
and dots refer to the numerical solution of the OBE for the 
probe and trigger group velocity, respectively. Parameters 
correspond to the perfectly symmetric situation considered 
in Fig. ini that is, F^^ = Fa = Ff ^ = F4 = F = 27r x 6 
MHz, ^1=^3^ 0.08F, ^2 = = F, 62=64 = 0, 
V i,j jij = = = 10"*F, and we have chosen 

M32 = M34 = 10"^" C -m, and N/V = 3.0 x lO" cm'^ Due 
to symmetry, one has perfect group velocity matching within 
a large interval of values for the detunings. 



equal probe and trigger group velocities in the case of the 
®^Rb five-level scheme specified above, at least within the 
entire EIT window. 

Finally, we note also that the problem with the pulse 
propagation, present in the asymmetric arrangement (c.f. 
Sec. Ill is absent in the symmetric arrangement. This 
is because the dispersion, and the related group velocity 
reduction have its origin in the linear part of the suscep- 
tibility for both, probe and trigger pulses. Hence there 
is no apparent singularity in the trigger group velocity, 
which in asymmetric case stemmed from the nonlinear 
origin of trigger dispersion. 



IV. CONCLUSIONS 

We have studied a five-level atomic system in two dif- 
ferent but related M-configurations. We focused on the 
nonlinear properties of the system and specifically on the 
conditions for the optimization of the cross-phase modu- 
lation between two weak fields of interest, which we have 
named probe and trigger fields. Both systems have been 
studied from a semiclassical point of view, i.e., by de- 
scribing all the fields in terms of their Rabi frequencies. 
We have seen that both linear and nonlinear properties 
are well described by an approach based on amplitude 
variables, which has been shown to reproduce well the 



scribed in the text, F^^ = F2 ~ 36 MHz, Ff ^ = F4 ~ 38 



H34 = 5.7 X 10" 
= 0.08F, 52 = S4 



C-m, 

0, V 



F^^' = = lO-^F, N/V = 3.0 X 10" cm-^ The 



FIG. 15: Group velocity of the probe and trigger pulses ver- 
sus the normalized detunings 5i/T = 83 /F. Full lines de- 
note the analytical predictions of Eqs. I53II and 15411 (the 
thick line refers to the probe and the thin line to the trig- 
ger). Circles and dots refer to the numerical solution of 
the OBE for the probe and trigger group velocity, respec- 
tively. Parameters correspond to the five-level scheme de- 
rived from the Di and D2 lines of the *^Rb spectrum de- 

^2 — 

MHz, fi32 = 1.27 X 10""^ C-m, 
Q4 = r, Q2 = 2.22F, fli = Q3 
ij lij 

asymmetry between the two dipole moment matrix elements 
has been compensated by adjusting the value of Q2- In this 
way group velocity matching is achieved within the entire EIT 
window. 



numerical solution of the exact optical Bloch equations 
describing the system. 

Both the asymmetric and the symmetric M scheme 
are able to provide a giant cross-Kerr modulation, which 
may be useful for many applications. Both M schemes 
can be seen as a "duplication" of the usual three-level A 
scheme at the basis of EIT, one for the probe and one for 
the trigger fields. In the asymmetric scheme, only the 
probe drives a significantly populated transition and a 
large cross-Kerr effect is obtained when either the probe 
or the trigger is slightly detuned from the two-photon 
resonance condition. The corresponding nonlinear phase 
shift, yielding for example the conditional phase shift of 
Eq. (P) of a quantum phase gate for photonic qubits, 
can become very large, especially when the probe and 
trigger group velocities, slowed down by EIT, become 
equal. In the asymmetric scheme, this group velocity 
matching can be achieved by properly adjusting the de- 
tuning and the intensity of the control field of the trigger 
A system. In the symmetric M scheme, the atomic pop- 
ulation is equally shared by the probe and trigger tran- 
sitions. Adjusting the corresponding parameters (Rabi 
frequencies, detunings) so that the two A systems be- 
come identical, probe and trigger experience the same 
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interaction with the atomic medium and group velocity two-photon resonance, so to remain still within the trans- 
matching is achieved automatically. In this case a signif- parency window. In fact, due to EIT, the susceptibility 
leant nonlinear cross-phase modulation is achieved only vanishes at all orders at the exact two-photon resonance 



if both A schemes are slightly and equally detuned from condition. 



APPENDIX A: OPTICAL BLOCK EQUATIONS - ASYMMETRIC CASE 

From Eqs. ijTJ) and H19|) . and using the definitions of Eqs. (|20f) - H22|l . one gets the following set of equations for the 
atomic populations an 

&ll = «rJi(72i - iriiCri2 + r4iO-44 + r2iO-22, (Ala) 
(722 = -i^l(J21 + i^lcri2 - 1^2(^23 + i^2'^32 " ^2(^22, (Alb) 
^33 = «f^30-43 - «f^3'^34 + «f^20-23 - «f^20'32 + r43Cr44 + r23Cr22, (Alc) 

(T44 = iQ^a:}4^ — iQ^a^s — ifl4^ai5 + iQ\a54^ — T4a44^, (Aid) 

(T55 = iri4(T45 — zri4Cr54 -I- r25Cr22 + r450-44, (Ale) 

and the following set of equations for the atomic coherences cry , i ^ j, 

P2 -|- ^12 

(712 = -i5i<7i2 + ifli{a22 - (Jii) - 1^2(^13 ^ o'i2, (A2a) 

0-13 = -«(5l20-13 + ifil(T23 - «^l3(Tl4 - if^2^i2 2"^^^^' (A2b) 

(7l4 = -«(5i3(Ti4 -I- zfii(T24 - «il30'l3 - »f^4f 15 - -^^^-^Y — -<7l4, (A2c) 

f^i5 = -^<^i4f^i5 + ^^1(^25 - i^lcri4 ^^15, (A2d) 

^ -k ^ -^2 H~ T23 

0-23 = i^2(723 + if^^CTig - ^17^(724 + iri2(CT33 - (722) ^ (723, (A2e) 

+ -^2 r4 -|- "^24 
6-24 = ^^230-24 — ^^3f^23 — ^^^4*^25 + ^^2*^34 + i^ldlA Cr24, (^^f) 

JJ'tj _[_ '~j/25 

0-25 = i62i<y25 + i^if^ib + i^2'^35 - i^Xcr24: ^"^ Cr25, (A2g) 

-1^4 H~ "^34 

(734 = -«<53(734 + iri3((744 - (733) + zrj2(724 - i^^40-35 ^ 1723, (A2h) 

^35 

0'35 = -«(5340-35 + ifi3CT45 + «^^20-25 - 1^4,0-34 2~^^^' 

(745 = Z(54(745 + irigCTss + iri4((755 - (744) - + 745 ^^^^ ^^2j) 

where we have also defined the composite detunings S23 = 62 — 63, S24 =62 — 63 + 64, and 634 = 63 — 64. 

APPENDIX B: OPTICAL BLOCH EQUATIONS - SYMMETRIC CASE 

From Eqs. 1)41) |l and (|19|) . and using the definitions of Eqs. I|2U() - (|22|I . one gets the following set of equations for the 
atomic populations an 

ail = «f^20-21 - «f^20'l2 + r4i(744 H- r2i(722, (Bla) 

(722 = «02(7i2 - 2ri20'21 - «f^l0'23 + «^^ia'32 - r2(722, (Bib) 

CT33 = «f^30-43 - «f^3^34 + «f^lO-23 - «^a-32 + r430-44 + r23(72, (Blc) 

(744 = lf23CT34 — 1573(743 — l5l4(745 -|- iri^a^4 — r4(744, (Bid) 

(755 = zr24(745 — ifi4(754 -I- r45(744 -I- r25(722, (Blc) 
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while the equations for the coherences are 



0'12 


1^2 ~t~ Tl2 

= -ld2<Jl2 + l^2{<722 - CTll) - lOiCTis fTi2, 

2 


(B2a) 


CTlS 


= i{6i - (52)ai3 + I^l2cr23 - i^lcTu - i^llai2 2"^^"^^' 


(B2b) 


(714 


= i{5i -62- S3)ai4 + iQ.\<j2A — 1^3(^13 — i^icrih — -^-^-^ctm, 


(B2c) 


O'lS 


= i{5i -82-63 + 5i)ai5 + iQ.\a25 - i^\crii ^"^is, 


(B2d) 


0'23 


= «^lCr23 + «f^20'l3 1^^3(724 + (0-33 Cr22) ^ ^23, 


(B2c) 


0'24 


= ^3)0'24 «fi30'23 l^i(J2b + t^2(^li + «f^l<''34 2 '^24, 


(B2f) 


0'25 


-|- "^25 

= + ^4 - '^3)0'25 + l^2<^\5 + ^^l'^35 " l^t<J2A 7, (^25, 

L 


(B2g) 


0'34 


= -«'530"34 + ^^^1(724 - «i^40'35 + l^3{aAi - (T33) ^-^--^Cr34, 


(B2h) 


O'SS 


= -«(<53 - 54)0'35 - l^^30'45 + «^^lCr25 - «f^4Cr34 '2'^^^' 


(B2i) 


0'45 


J , \ , r>* r4 + 745 
= l()40-45 + 11^4(0-55 (744j + 14^30-35 (745- 


(B2j) 
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